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Abstract 

We realize, by Monte Carlo event generator methods, the exact 0(a) YFS exponenti- 
ated calculation of e + e~ — > W + W~ '(— > fif[ + fif£) at and beyond LEP2 energies, where 
the left-handed parts of fa and // are the respective upper and lower components of an 
SU2L doublet, i = 1,2. Our calculation is gauge invariant from the standpoint of its 
radiative effects and the respective YFS Monte Carlo event generator YFSWW3, wherein 
both Standard Model and anomalous triple gauge boson couplings are allowed, generates 
72(7) radiation both from the initial state and from the final W + W~ . Sample Monte Carlo 
data are illustrated. 
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1 Introduction 



The processes e + e~ — > W + W~ +^(7) — > 4fermions + n('y) at and beyond LEP2 energies 
are of considerable interest for the LEP2 and NLC physics programs. When calculated 
and measured to sufficient precision, they provide windows to several important avenues 
for verification and tests of the SU2L x U\ model of Glashow, Salam and Weinberg |J of 
the electroweak interaction. In Ref. || for example, the goal for the theoretical precision 
of the understanding of these processes was set at 0.5% for the LEP2 physics program. In 
the following discussion, we will present calculations of these processes which feature the 
most precise electroweak radiative corrections which have been published to date, that 
of an exact 0(a) YFS || exponentiated Monte Carlo event generator in which multiple 
photon radiatve effects from both initial and W + W~ states are realized on an event- 
by-event basis as infrared singularities are canceled to all orders in a. We denote the 
respective event generator by YFSWW3. We will comment below on its physical and 
technical precision relative to the 0.5% precision goal of Ref. 0. 

A few comments are in order regarding the relationship between the calculations pre- 
sented below and those presented in Ref. ||]. In Ref. ||, we presented the second order 
leading-log YFS exponentiated calculation of the processes under discussion here and in- 
cluded for the first time the complete YFS formfactor effect for the radiation from the 
W + W~ themselves, in a gauge invariant way || §]. With the advertised || precision 
goal of 0.5% for the final LEP2 theoretical precision tag on the signal prediction for the 
W + W~ pair production, it is important to treat the complete 0(a) exactly in the pres- 
ence of our YFS exponentiation, including the exact pure weak corrections at 0(a), as 
these may enter at the level of this goal. It is for this reason that we are motivated to 
carry-out our work in this paper. 

Specifically, what we will do is the following. In the next section, we shall describe the 
extension the analysis in Ref. @] to include final state multiple photon radiation at the 
leading-log (LL) level. This will lead to the development of the new YFS exponentiated 
Monte Carlo event generator YFSWW3, involving as it will the extension of the YFS3 
event generator in Ref. 0] for the process e + e~ — > f f + 71(7) to the process e + e~ — > 
W + W~ + 71(7), where the W-pair is then allowed to decay. This will set the stage for 
our work in section 3, wherein we will include in the YFSWW3 event generator the exact 
0(a) results for the hard photon residuals i — 0, 1 in the language of Ref. (this 
language will be reviewed briefly below) basing ourselves on the exact virtual electroweak 
0(a) results and the exact O(o) bremsstrahlung results of Refs. ||. In section 4, we will 
then present some sample Monte Carlo data to illustrate the size of the various levels of 
approximation that we access in our work, with some focus on what is currently available 
in the literature and how our results relate to it. Section 5 contains some concluding 
remarks. 
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2 Multiple Initial + Final State Photon Radiative 
Effects 



In this section we discuss the relevant aspects of our YFS Monte Carlo methods as 
they pertain to the problem of extending our YFS3 MC in Ref. for the processes 
e + e~ — > ff + 72(7) i n which both initial and final state YFS exponentiated multiple pho- 
ton radiation is simulated to the W + W~ pair production processes of interest to us here. 
We carry-out the respective extension as well to arrive at a new MC event generator, YF- 
SWW3, in which not only is the full YFS formfactor effect with the VT-pair contribution 
included as in YFSWW2 || but also real soft final state 71(7) radiation from the W-pair 
itself is calculated. Such results have not appeared elsewhere. On the way we also set our 
notation and define our kinematics for the work in the remaining sections. 

Specifically, the processes of interest to us here are illustrated in Fig. [TJ, where we have 
given our kinematics. 




Figure 1: The process e + e — > W + W Afermions = fx + f[ + f 2 + f' 2 , where \jtj > « = 1, 2, 

are SU2L doublets. Here, pa is the 4-momentum of A, A = fi,f-, Pi(qi) and ^2(^2) are the 4- 
momenta of e + (e~) and Vt /+ (V4 / ~) respectively. We use the notation Cl = Pl C = |(1 — 75) C for 
any C. 

We consider e + + e~ — > W + + W~ + 71(7) — > Afermions + 71(7) at CMS energies yfs > 
2Mw In Ref. 0, we have shown that our YFS exponentiated Monte Carlo algorithms in 
Refs. H for the processes e + e~ — > / + / +71(7), where / is a fundamental fermion in the 
SU 2 l x U\ theory [|IJ, has a gauge invariant extension to e + + e~ — > VT" 1 " + W 7- + 71(7) — > 
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Afermions + 71(7). In Ref. [Q we presented the extension for our YFS2 initial state 
multiple photon Monte Carlo as the YFSWW2 multiple photon Monte Carlo in which 
the full YFS form factor was taken into account. Here, we present the corresponding 
extension of our YFS3 initial + final state multiple photon Monte Carlo in which mutiple 
photon radiation from the Ws themselves is included in the in the lowest level of the 
Monte Carlo algorithm. We refer to this extension of YFS3 to the H^-pair production 
process as YFSWW3 and it is available from the authors . 



Specifically, on applying the results in Ref. [§J for extending, to IV-pair production in 
a gauge invariant way, our YFS formula for the process for e + e~ —>/ + / + ^(7) to our 
formula for the cross section in Ref. JtJ, we arrive at the gauge invariant YFS formula 
for the respective VT-pair production cross section in which multiple initial state and 
multiple final state radiation is calculated, where we use the ratio of T w /M w to justify 
treating the W^-pair state as a "final" state on their mass shell insofar as these radiative 
effects are concerned. We then realize this formula by Monte Carlo methods in complete 
analogy with what is done in Refs. [|, |7|, using as our Born cross section that given 
by Ref. [ll]. The resulting Monte Carlo, YFSWW3, at high energies compared to My/ 



features final state leading-log 0(a 2 ) radiative effects as well as the initial state leading-log 
0(a 2 ) radiative effects in YFSWW2. Thus, at NLC energies, these final state radiative 
effects may be important as we shall illustrate presently. In addition, the YFSWW3 final 
state 71(7) radiation allows us to check the technical precision of the /?o-level 71(7) soft 
radiation in YFSWW2 since the final state 71(7) radiation is already present in the low 
level background Monte Carlo in YFSWW3 whereas in YFSWW2 it is not. We will also 
illustrate this presently. 



Ecu [GeV] 


ISR 


+ Coul. corr. 


+ Y'-corr. 


± WW-md. 


175 


0.4906 ± 0.0002 
0.4898 ± 0.0002 


0.5046 ±0.0002 
0.5037 ±0.0002 


0.5053 ±0.0002 
0.5048 ±0.0002 


0.5053 ±0.0002 
0.5048 ± 0.0002 


190 


0.6060 ±0.0007 
0.6034 ±0.0007 


0.6193 ±0.0007 
0.6166 ±0.0007 


0.6217 ±0.0007 
0.6195 ±0.0007 


0.6219 ±0.0009 
0.6197 ±0.0009 


205 


0.6359 ± 0.0008 
0.6315 ±0.0008 


0.6480 ±0.0008 
0.6436 ±0.0008 


0.6514 ±0.0008 
0.6476 ±0.0008 


0.6516 ±0.0010 
0.6475 ±0.0010 


500 


0.2910 ±0.0003 
0.3538 ± 0.0004 


0.2946 ±0.0003 
0.3582 ±0.0004 


0.2970 ±0.0004 
0.3591 ±0.0004 


0.2954 ±0.0004 
0.3571 ±0.0005 



Table 1: The results of the 10 5 (except for Ecu = 175 GeY, where it is 10 6 ) statistics sample 
(unweighted events) from YFSWW3 for the total cross section a \pb]. The upper results at each value 
of energy are for the Standard Model couplings constants, while the lower ones are for anomalous 
couplings constants (5k = 5X = 0.1). See the text for more details. 

More precisely, in Table [l], we show the result of our YFSWW3.0 simulation of the 
process in Fig. [I] with the cs ± ev e final state, for CMS energies 175 GeV, 190 GeV, 205 
GeV and 500 GeV which have LEP2 and NLC in mind. We present results according 
to the notation in Ref. ||] so that ISR denotes initial state 71(7) radiation, "Coul. corr." 
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denotes that the Coulomb correction after the fashion of Ref. is included, "Y'-corr." 
denotes that the full YFS form-factor effect already featured in Ref. [[| is turned on, and 
finally 'WVT-rad." denotes the new feature of YFSWW3.0 in which the real soft 7^(7) 
bremsstrahlung from the W-pair is included in the simulation (already at the level of 
the low level Monte Carlo algorithm in the sense of Ref. || for example). In addition, 
the upper and lower results in each entry in the table correspond to the case of Standard 
Model and anomalous (5n = 5X = 0.1) WWV couplings in the notation of Ref. [[11]]. Thus, 
the comparison of the last two columns shows that the effects of the final state ^(7), while 
negligible at LEP2 energies, is significant at NLC energies. Comparison of the first three 
columns in the table with the analogous results in the respective Table 1 of Ref. shows 
that indeed we have a very good agreement between YFSWW2 and YFSWW3.0 in the 
calculation of effects in which the real final state radiation should not be important. This 
represents a very good technical precision check on the two calculations. Moreover, it 
shows that YFSWW3.0 is an excellent starting point for developing the exact 0(a) YFS 
exponentiated Pi-level Monte Carlo event generator calculation in which one has an exact 
0(a) calculation in the presence of initial + final state 0(a 2 ) leading-log radiation with 
YFS exponentiation - such a calculation has not appeared elsewhere. To this we now 
turn in the next section. 



3 Exact 0(a) YFS Exponentiated WW-Pair Produc- 
tion: YFSWW3.1 

In this section, we develop the exact 0(a) realization of the hard photon residuals j3 n , n = 
0, 1 in our YFSWW3 Monte Carlo event generator, where we refer to Refs. £J |J for a 
precise definition of these residuals. We start with /3 . 

For the construction of the exact 0(a) electroweak result for /3q, which we denote by 
Pq , we first note that 

-pM = da one - loop /dQ - 2aMBda B /dtt (1) 

where the YFS virtual infrared function B for the process e + + e~ — > W + + W~ is defined 
in Ref. 01], where the cross section das/dQ is the respective Born cross section, and where 
da one ~ loop / dQ is the respective exact one-loop correction to the cross section given by the 
results of Ref. |£| . We have implemented the formula in (J!]) in our YFSWW3 Monte Carlo 
event generator to arrive at a realization of the exact 0(a) electroweak result for the hard 
photon residual j3 . 

Turning now to the exact 0(a) result for the hard photon residual which we denote 
by Pi i we first note that 

~p[V = da B1 /kdkdQ 1 dQ - S(k)da B /dtt (2) 
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where S(k) is the respective YFS real infrared function defined in Ref. Q and da B1 / kdkdQ^dQ 
is the respective exact 0(a) bremsstrahlung cross section^] from Ref. || for the photon 
into the phase space element kdkdQ^ when the W~ is produced into the solid angle dfl, 
for example. We have implemented the result (0) into our YFSWW3 Monte Carlo as well 
and the resulting version of it, including both the results in ([!]) and in (fj) is version 3.1: 
YFSWW3.1. In the next section, we illustrate some of its applications. 

4 Illustrative Results from YFSWW3.1 

In this section we illustrate the application of the exact 0(a) YFS exponentiated Monte 
Carlo event generator YFSWW3.1 to the WW-pair production and decay at LEP2 and 
NLC energies. We continue to use the cs + ev e 4-fermion final state for definiteness. 

Specifically, we have recorded in Table a summary of the sizes of the various approx- 
imations which we have realized in YFSWW3.1 including the exact 0(a) YFS exponen- 
tiated result. 



Ecu [GeV] 


a [pb] 


(a p 2 rag -a Q )/a 




{of - a^/ao 


(°f - °f )/*o 


161 


0.1768 


-24.6% 


+0.11% 


-0.80% 


+0.21% 


175 


0.5891 


-15.4% 


+0.12% 


-1.31% 


-0.007% 


190 


0.6792 


-10.0% 


+0.11% 


-1.70% 


-0.22% 


205 


0.6850 


-7.1% 


+0.10% 


-2.20% 


-0.61% 


500 


0.2710 


+4.4% 
+4.4% 


+0.19% 
+0.19% 


-4.65% 
-5.05% 


-3.08% 
-3.08% 



Table 2: Various contributions (approximations) to the YFS exponentiated WW cross section as 
a fraction [in %] of the Born cross section a?, ra9 denotes the so-called pragmatic 0(a 2 ) cross 
section (see the text for more details), a\ L and a 2 L are the LL approximations of 0(a 1 ) and 0(a 2 ), 
respectively, af x is the exact 0(a 1 ) result, while a°y is the so-called Improved Born Approximations 
(see the text for more details). The results in the lower line for Ecu = 500 GeV correspond to the 
situation when LL QED corrections for the W + W~ state are also included in the corresponding LL 
contributions. 

We show results for CMS energies E CM of 161, 175, 190, 205 and 500 GeV. The 
definitions of the YFS exponentiated cross sections in the table are as follows: a^ ag 
denotes that the cross section contains the exact O(o) correction and the 0(a 2 ) initial 
state LL correction at the level of the W-pair production process; a^ L denotes that the 
cross section contains initial state LL QED radiative corrections through 0(a n ); o e * 
denotes that the cross section contains the exact radiative correction through 0(a n ); 
a°f indicates that the cross section contains an approximate treatment of the radiative 
corrections through 0(a n ) (so-called Improved Born Approximation) and we only discuss 
the case n — 1 where we have implemented the approximate cross section of Ref. |13j . The 

1 We would like to thank K. Kolodziej for providing us with routines for the 0(a) hard bremsstrahlung 
matrix element. 
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second entry for Ecm = 500 GeV corresponds in each column to a different treatment of 
the LL QED radiation in which both initial and final state LL QED effects are included 
- it gives us an estimate of the size of the sub-leading QED radiative effects at the higher 
orders when it is compared with the first entry for example. The cross sections shown in 
the table correspond to simulations of 10 6 (weighted) events at all energies except for the 
Ecm = 500 GeV results, where 10 7 events were simulated. What we see in the table is that 
the 0(a 2 ) LL effects are well below the 0.5% targeted precision tag of the Physics at LEP2 
"WW Cross-Sections and Distributions" Physics Working Group || for all LEP2 energies. 
For the 500 GeV case, the high precision objectives of certain NLC physics issues would 
clearly necessitate inclusion of these second order LL effects. The difference between the 
exact 0(a) result erf and the first order LL result a\ L shows that even at LEP2 energies, 
the non- leading QED and pure EW corrections are important Jl3j; and, as advertised in 
Ref. [T3|, the use of the approximation cr" p at LEP2 energies would give results within 



the desired 0.5% precision tag except perhaps at the highest LEP2 CMS energy of 205 
GeV, where the difference between the exact result and the a^ p is 0.61%, just beyond the 
desired tag. Thus, for all but the highest LEP2 energies, the cr" p could be used to realize 
a gain in CPU time for event generation of a factor of ~ 30 without an unacceptable loss 
in precision relative a 0.5% precision tag. At the NLC type energy Ecm = 500 GeV the 
approximation <r" p is not supposed to be accurate and our last column bears this out, 
although we see that it is better than the o\ L at such high energies. We are encouraged 
that our best result, <rf ra9 , is not very sensitive to the inclusion of the final state LL QED 
(a 2 ) effects at 500 GeV; this is not true of crf i , as expected. We stress that, in the final 
state, L = ln(s/Myy) — 1 is only 2.67 for Ecm — 500 GeV so that the LL series in the final 
state is one of several effects of comparable significance in comparison to the LL initial 



state series at 
effects. 



2 ) and we show the last row in the table to illustrate the size of such 



Ecm [GeV] 




(# 1J -#!L)/*o 


161 


-0.85% 


+0.04% 


175 


-1.24% 


-0.07% 


190 


-1.57% 


-0.13% 


205 


-2.06% 


-0.15% 


500 


-5.38% 
-6.05% 


+0.73% 
+1.00% 



Table 3: Differences between exact and LL approximated non-infrared contributions to the YFS 
exponentiated WW cross section as a fraction [in %] of the Born cross section gq. The results in 
the lower line for Ecm = 500 GeV correspond to the situation when LL QED corrections for the 
W + W~ state are also included in the corresponding LL contributions. 



Since we have incorporated the exact 0(a) correction in our calculation, it is instruc- 
tive to look at the size of its effect in comparison to the LL QED radiative correction 
insofar as our YFS hard photon residuals are concerned. We show this in Table |3|. In 
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the table, the notation is such that is the 0(a m ) result for the hard YFS photon 
residual /3 n and the subscript LL denotes that it is computed in the LL approximation 
(in the initial state LL approximation to be precise). What we see in the table is that 
the main effect of the exact result in comparison to the LL one is realized in the hard 
photon residual /?o; the effect on (3\ in comparison to the LL approximation is well below 
the 0.5% precision tag of Ref. |2| for LEP2 energies. For the NLC energy 500 GeV, we 
see that, while the main difference between the exact result and the LL approximation is 
still in (3 , that in j3\ is also significant for high precision NLC physics objectives. 

These results on the total cross sections are then analyzed in more detail in Figs. 0, 
H, f|, wherein we show the W~ angular distributions corresponding to the respective Eqm 
values in Tables 1-3. For each value of Ecu we show the W~ angular distribution at the 
Born level and at the O^a 2 )?™ 9 level (crf raff level) as well as the differences Ada{9) 2 -i = 
da% L (9)-daf L (9), Ada(9) ex ^ = daf x (9)-daf L (9), and Ada(9) ex _ ap = daf '(9) - dtr? \9) 
in ratio to the respective Born level differential cross section dao(9), where Fig. 1 con- 
tains this information for Eqm = 161,175 GeV, Fig. 2 contains it for Eqm = 190,205 
GeV, and Fig. 3 contains it for E CM = 500 GeV. In the figures, we denote by 0(a n )^ xp 
the YFS exponentiated cross section a^ L level, by 0(a n ) exp the YFS exponentiated cross 
section level, and by 0(a n )^ rox the YFS exponentiated cross section a°f level. The 
differential cross section results are fully consistent with the total cross section results 
shown above. Moreover, they show us that at all energies, the radiative corrections are 
significant throughout the entire angular distribution. At LEP2 energies, the angular de- 
pendence of the ratio Ada (6 , ) 2 -i /da (9) is essentially flat; at the NLC energy 500 GeV, it 
is nontrival and strongly varying in the backward direction. The ratios Ada{9) ex -\ / dcr (9) 
and Ada(9) ex _ ap /dao(9) are smoothly varying in 9 at LEP2 and at NLC energies; at all 
energies, in the region where the cross sections are largest, the latter is smaller in mag- 
nitude than the former, as expected; away from this region, especially near the backward 
direction, the opposite is true. Of course, at the NLC energy 500 GeV, one does not 
expect the cross section a\ p be as accurate as it is at the LEP2 energies for which it was 
optimized and this is borne-out both in the 9 dependence of the ratio Adcr(9) ex _ ap /d<To(9) 
as well as in the total cross section results in Table 2. Not shown in the Figs. 2-4 are the 
corresponding results for the initial + final state LL YFS exponentiated approximation 
at 500 GeV; its curves are entirely similar to those shown in Fig. 4 so that we do not 
present them separately here. 

Our conlusion is that the exact results show that for the targeted precision of the 
Physics at LEP2 Workshop the LL approximation is inadequate, both for the total cross 
section and for the differential cross section ; the non-leading corrections must be taken 
into account. In YFSWW3.1, we have the only multiple photon YFS exponentiated 
amplitude based calculation of these effects in which the respective infrared singularities 
are canceled to all orders in a and in which the corrections themselves are realized on an 
event-by-event basis. 
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Figure 2: The W~ angular distributions for two LEP2 energies: yfs = 161 GeV (left picture) and 
y/s = 175 GeV (right picture). The upper parts of the pictures show the differential cross sections 
in Born (thin-lines) and in 0(a 2 )exp 9 (thick-lines) approximations, where the latter one denotes the 
YFS exponentiated cross section including the exact 0(a 1 ) matrix element together with the 2-nd 
order LL ISR terms (see the text for more details). In the lower part of each picture the following 
contributions to the YFS exponentiated cross section (divided by the corresponding Born differential 
cross section dao) are presented: the 2-nd order LL ISR contribution (oo), the 1-st order non-LL 
contribution including electroweak radiative corrections (oo) and the difference between the exact 
0(a l ) matrix elements and the so-called Improved Born Approximation (**). The last one shows 
the quality of IBA for the differential cross section as a function of the VF-boson polar angle. 

5 Conclusions 



In this paper we have presented and illustrated the first ever 4-fermion Monte Carlo event 
generator for W + W~ pair production and decay in which the full EW 0(a) correction 
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Figure 3: The same as in Fig. || but for other LEP2 energies: y/s = 190 GeV (left picture) and 
V's = 205 GeV (right picture). 



is calculated in the presence of amplitude based YFS exponentiation in which infrared 
singularities are canceled to all orders in a. The resulting Monte Carlo event generator, 



YFSWW3.1, is available from the authors [1C 



As we have illustrated by showing the results of simulations at several LEP2 ener- 
gies and at the NLC energy 500 GeV, the nonleading 0(a) EW corrections are indeed 
important and must be taken into account, both for the overall normalization and for 
the detailed angular distributions. At LEP2 energies, we showed that the approximate 



representation in Ref. [[T3| of these non-leading corrections does indeed stay with 0.61 of 
the exact 0(a) EW result in the presence of our YFS exponentiation as we realize it by 
Monte Carlo methods; for the NLC energy, this is no longer the expected. We 

stress that these last remarks apply both to the total cross sections and to the differential 
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Figure 4: The same as in Fig. ||| but for the NLC energy: ^/s = 500 GeV. 



cross sections in the regions where they have their largest values. Our j3\ level second 
order LL YFS exponentiated results are thus a significant improvement over the (3q level 
results in Ref. 0]; it is also a significant improvement over the other available Monte Carlo 
event generators |T4] , |T5| in the literature; for, none of them feature a 4-fermion, amplitude 
based, exact EW 0(a) YFS exponentiated calculation in which infrared singularities are 
canceled to all orders in a, on an event-by-event basis as we do in YFSWW3.1. We look 
forward with excitement to the many applications of our work. 
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